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Abstract 

Many observers in the past gave detailed descriptions of the telescopic as- 
pect of Venus during its extremely rare transits across the Solar disk. In 
particular, at the ingress and egress, the portion of the planet's disk outside 
the Solar photosphere has been repeatedly perceived as outlined by a thin, 
bright arc (" aureole" ) . Those historical visual observations allowed inferring 
the existence of Venus' atmosphere, the bright arc being correctly ascribed to 
the refraction of light by the outer layers of a dense atmosphere. On June S*'^, 
2004, fast photometry based on electronic imaging devices allowed the first 
quantitative analysis of the phenomenon. Several observers used a variety 
of acquisition systems to image the event - ranging from amateur-sized to 
professional telescopes and cameras - thus collecting for the first time a large 
amount of quantitative information on this atmospheric phenomenon. In this 
paper, after reviewing some elements brought by the historical records, we 
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give a detailed report of the ground based observations of the 2004 transit. 
Besides confirming the historical descriptions, we perform the first photomet- 
ric analysis of the aureole using various acquisition systems. The spatially 
resolved data provide measurements of the aureole flux as a function of the 
planetocentric latitude along the limb. A new differential refraction model 
of solar disk through the upper atmosphere allows us to relate the variable 
photometry to the latitudinal dependency of scale-height with temperature 
in the South polar region, as well as the latitudinal variation of the cloud-top 
layer altitude. We compare our measurements to recent analysis of the Venus 
Express VIRTIS-M, VMC and SPICAV/SOIR thermal field and aerosol dis- 
tribution. Our results can be used a starting point for new, more optimized 
experiments during the 2012 transit event. 

Keywords: Venus, Venus atmosphere. Planet transits 



1. Introduction 



Since the mid 18*'^ Century, observers have reported unusual features of 
the telescopic image of Venus near the inferior conjunction, promptly at- 
tributed to its atmosphere. Some of them are at the reach of modest in- 
struments (although at small angular distance from the Sun), such as the 
cusp extension first described by Schroeter ( 1791[ ), which tends to transform 
the thin crescent of Venus into a ring of light (e.g. Russell, 1899 DoUfus & 



Maurice 1965). 



One of the most relevant features pertaining to ground-based studies of 
the Venus atmospheric structure has been observed only during transits, 
close to the phases of ingress (between 1st and 2nd contact) ot egress (be- 
tween 3rd and 4th contact) as a bright arc outlining - in part or entirely - 
the portion of Venus' disk projected outside the solar photosphere. Tradi- 
tionally the first account of this phenomenon was attributed to Mikhail V. 
Lomonosov (1711-1765) who reported his observations of the transit at St. 
Petersburg Observatory on May 26, 1761 (Marov, 2005). Actually the poor 



performance of his small refractor hints that most probably other observers 
(such as Chappe d'Auteroche, Bergman, and Wargentin) were the first gen- 



uine witness on the same date (Link, 1969 Pasachoff and Sheehan, 2012). 
However, Lomonosov correctly attributed the putative phenomenon to the 
presence of an atmosphere around the planet, refracting the sunlight in the 
observer's direction. 
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In the following, adopting a denomination widely used in the historical 
accounts, we will often call this arc "aureole". Since both the aureole and 
the cusp extension occur close to the planet terminator, they are also col- 
lectively known as "twilight phenomena". For detailed historical reviews the 
interested reader can refer to Link (1969) and Edson (1963). Of course, in 
this context we neglect the initial motivation for the observation of transits: 
the determination of the Astronomical Unit by the solar parallax (proposed 
by E. Halley in 1716), whose interest is purely historical today. 

Venus transits are rare, as they occur in pairs 8 years apart, each pair 
separated by 121.5 or 105.5 years, alternating between descending node 
(June pairs: 1761/1769, 2004/2012), and ascending node (December pairs: 
1631/1639, 1874/1882, 2117/2125). As a consequence, data concerning the 
aureole are correspondingly sparse and, up to the last event, they have been 
obtained by simple visual inspection, mainly by using refracting telescopes of 
modest aperture (typically up to 15-20 cm). This limitation was mainly due 
to the constraint of organizing complex expeditions including the delicate 
transportation of all the instruments. 

The 2004 event represents a giant leap in the observation of Venus tran- 
sits, as the modern imaging technologies available allow for the first time a 
quantitative analysis of the atmospheric phenomena associated to the transit 
of Venus. 

We performed measurements of the aureole on the original images ob- 
tained through several different instruments, and compared them both with 
a simple refraction model and with observations obtained in the past. This 
work summarizes the aspect of Venus close to the Sun's limb during the 
June 8, 2004 transit as observed by ground-based instruments. A recent pa- 
per (Pasachoff et al. , 2011) deals with imaging using NASA's then operating 
Transition Region and Coronal Explorer solar observatory (TRACE). 

Since the ground-based observations were not specifically organized be- 
forehand, data at our disposal are rather heterogenous. In order to bridge the 
gap between past visual observations through small telescopes and today's 
technologies, we decided to consider accounts obtained both with profes- 
sional instruments and low-cost amateur telescopes, either by CCD imaging 
or by direct image inspection by experienced observers. In fact, as shown 
in recent studies of distant solar system objects based on stellar occultation 
campaigns (e.g. Widemann et al., 2009) while CCD imaging offers today the 
most valuable quantitative measurements, small telescopes and visual obser- 
vations allow a significant increase of constrains on the phenomena - and in 
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the case of Venus, the most direct comparison to past reports. The results 
obtained from the analysis of the most significant image sets, representing a 
certain variety of instrument size and quality, are illustrated in this paper. 

The paper is organized as follows. First, we describe the conditions of 
the 2004 event, and the reconstruction of limb geometry (section |2]). We 
then provide an extensive review of the measured spatial and temporal vari- 
ations for the brightness of the aureole in wide or narrow-band photometry 
(section [3] and |4]). We then address the basic physical principles of the at- 
mospheric differential refraction model producing the aureole (section [5| and 
we use it for the interpretation of the observations ([6]). In Section [7| the 
modeling of the aureole is compared to recent analysis of the Venus Express 
observations (VIRTIS-M, VeRa and SPICAV/SOIR) regarding the thermal 
field and cloud-deck altitude and haze distribution pertaining to this study, 
as well as ground-based mid-infrared spectroscopy of non-LTE CO2 emission. 
The comparison is discussed. 



2. Geometry of the transit in 2004 



For simplicity and following Link (1969) hereinafter we call "phase" (/) 



of the event the fraction of Venus diameter external to the solar photosphere. 
A value f=0 corresponds to the planet entirely projected on the Sun, tangent 
to its limb. When /=0.5, the planet centre will be exactly on the solar limb, 
and so on. 

In Fig. [1] the orientation of the disk of Venus relatively to the solar limb 
is given, both for ingress and egress phases. In both cases the temperate 
latitudes are tangent to the Sun for / = and / = 1. At first contact the 
South pole of Venus remains projected longer on the sky, while the North 
pole is the first one to enter on the solar disl^ The sequence is inverted 
between third and fourth contact, such that it is always the South pole to 
be observed externally to the Solar disk for a longer time. 

The total limb crossing for the disk of Venus lasted 18.9 minutes and the 
apparent radius of the planet was 28.9 arcsec. 



^In the following we will always use the lAU convention i.e. the North pole is the one 
lying on the northern side of the ecliptic. 
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Figure 1: Sketch of Venus' disk orientation at the ingress (left panel) or egress (right) of the 
transit. The greyed area corresponds to the Sun's photosphere at the epoch when the Venus 
disk is externally tangent to the Sun (first and fourth contact: ti and t4 respectively). 
The Solar limb is also indicated at the second and third contact (labeled t2 and t^) and 
when projected on the center of Venus. The solid arrows indicate the direction of the 
apparent motion of the planet relative to the Sun. The vertical dashed line corresponds 
to the sky North-South direction, while the dash-dotted one represents the sky-projected 
rotation axis of Venus. 
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3. Observations and measurements 



The European observers providing tlie data sets described further on had 
particularly favourable conditions around the end of the event, while the 
ingress of the planet on the solar disk was observed at much higher airmass, 
i.e. at low elevation above the horizon. However, visual observers under 
good sky conditions and employing a magnification higher than ~150x had 
no particular difficulty in identifying the bright aureole outlining the Venus 
disk between 1st and 2nd contact, while it was crossing the Solar limb (i.e. 
for f<l). Skilled observers immediately noticed the non-uniform brightness 
of the aureole along the planet disk (a high-quality drawing by an expert 
amateur observer is shown in Fig. [2| . 




Figure 2: Drawings of Venus at the beginning (three leftmost panels) and end of the event 
(remaining two panels), as seen visually by an expert amateur astromer (Mario Frassati, 
Italy) through a 20 cm Schmidt-Cassegrain telescope. Courtesy of the archive of Sezione 
Pianeti, Unione Astrofili Italiani. 



For a quantitative analysis of this phenomenon we rely on images ob- 
tained by CCD cameras through different telescopes. Given the casual na- 
ture of most images of the aureole that have been produced during the event 
through instruments of all sizes, not all of them are suited for a compre- 
hensive analysis. We thus selected the representative sample of observations 
presented in the following, with appropriate image quality and information 
content. Table [T] summarizes the actual contributions that were collected 
and analyzed, as presented in the following. 
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Table 1: Measured observations. The availability of ingress ("In" column) or egress ("Ex") 
images is given. The number of measured images follows. These are final images, as a 
result of an averaging process when needed. The final column identifies the filter name 
and/or their central wavelength/width. See case by case detailed explanations in the text. 



Observer 


Instrument 


Site 


In 


Ex 


N. images 


A/FWHM (nm) 


J. Arnaud 


Themis 


Tcncrife (Spain) 




X 


50 


V 


L. ComoUi 


20 cm Schmidt-Cass. 


Tradate, Italy 


X 


X 


50 


V 


P. Suetterlin 


DOT 


La Palma (Spain) 




X 


21 


430.5/1.0, G 










X 


21 


431.9/0.6, Blue cont. 










X 


21 


655.0/0.6, Red cont. 










X 


21 


396.8/0.1, Ca II H 


S. Rondi 


50 cm TourcUc rcfr. 


Pic du Midi (France) 


X 




32 


NaDl/10 



Images have undergone a standard calibration process (dark current and 
flat-field correction). The brightness of the aureole, when present, was mea- 
sured. Fig. [3] shows an example of CCD image sequence between 3rd and 
4th contact on June 8, 2004. The bright arc or aureole is clearly visible in 
the three frames obtained at an increasing distance of the Venus disk center 
from the Sun. For illustration purposes a contrast enhancement is applied. 
The inhomogeneity of the aureole is clearly visible on the motionward limb 
of Venus. 




Figure 3: Three images concerning the final phases of the transit obtained by the DOT 
telescope in the "G" band. Each image has been altered in contrast (gamma — 0.4) in 
order to better show the thin and faint arc of light outlining (partially or completely) the 
Venus limb. The orientation is the same as in Fig. [l] The bright saturated area is the 
solar photosphere. The progressive reduction in extension and brightness of the aureole 
with time is clearly visible. 

The arc's photometric profile was obtained by the integration of the signal 
in a ring centered on the planet, outlined in Fig. |4} The full ring is divided 
into sectors of identical angular extension (as seen from the center of the 
Venus disk), each corresponding to a flux integration area. Among them, 
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only those sectors projected on the sky are then considered. The background 





Figure 4: Scheme of the procedure adopted for the measurements and implemented in 
a custom software. Here black is used to represent the sky background and the disk 
of Venus. The four white circles define three concentric rings, with the aureole entirely 
contained in the central one. The nearby external and internal rings are used for evaluating 
the background. Each ring is divided in sectors of equal angular size. One of them 
is represented for the aureole (solid rectangle) and the corresponding background areas 
(dotted). The dashed angle represents the position angle (relative to the celestial North) 
used for representing the brightness profiles in following plots. The dashed-dotted line 
represents the orientation of the rotation axis of Venus. 



is evaluated on two rings, one on the inside and the other on the outside of 
the main ring. Their outer and inner edges (respectively) are also plotted 
in Fig. |4j Each background ring is divided in sectors as the measurement 
annulus. This way, each patch on the main annulus is accompanied by two 
background patches at the same position angle (relative to the Venus disk 
center) but on opposite sides. Their average flux, weighted by their surfaces, 
yields the final background value to be used for the given measurement patch. 
In the images we selected, the solar disk is not saturated, so a normal- 
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ization relative to a selected photosplieric region is possible. We decided to 
use the value of brightness measured at 1 Venus radius from the solar limb. 
On the images, the geometric limb is estimated by the analysis of the radial 
brightness profile of the solar disk. The photometric normalization value is 
provided by the average of the flux in a square window of the same surface as 
the sectors used for measuring the aureole brightness. A final normalization 
is done for deriving the plotted values, corresponding to the flux coming from 
an unresolved, thin segment of aureole 1 arcsec in length, normalized to the 
brightness of a square patch of photosphere having a surface of 1 arcsec^. 

In the following, we present technical details of each data set, with the 
corresponding aureole luminosity derived. The error bars given in the plots 
are computed by considering the Poissonian statistics of the photometric 
noise. We indicate by Ng and Ni, the number of photoelectrons coming 
from the aureole and the background, respectively, in each arc bin. The 
uncertainty that we consider is a = {Ng + Nb)^-^. In our case, Ng, thus 

the only relevant contribution to noise comes from the background, and is a 
function of the distance of the arc element from the solar limb. The camera 
read-out-noise is cleary smaller than the background contribution and can be 
discarded as well. The same applies to noise associated to the dark current, 
since exposure times are extremely short. 

From a practical point of view, we have verified that a appears to be fairly 
constant over the entire arc, growing of about 50% only for measurement 
points very close to the solar limb. For this reason, we choose to show only one 
error bar, common to all curves, except in the case of observations obtained 
at the DOT telescope which are used for modelling the aureole flux. In the 
related plots, two error bars separately for low and high / values, provide an 
additional information on the amount of noise variation in time. 



3.1. Tourelle Solar Telescope - Pic du Midi, France 

Only the ingress has been imaged by this station, since clouds have pre- 
vented a clear view of the final phases. Given the image scale, the radius of 
Venus is 142.5 pixels. For a set of selected images in which the aureole is 
well visible at the visual inspection, its integrated flux was determined us- 
ing different widths of the measuring ring, from 2 to 14 pixel. The resulting 
"growth curve" (Howell, 1989) shows the convergence toward the inclusion of 
the complete flux coming from the arc. The flux was considered to be com- 
pletely inside the ring at a width of 9 pixel, then used for all measurements. 
The arc was sampled at angular steps of 9 degrees. 
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Also, it was found that due to the faintness of the aureole, a useful im- 
provement of the SNR was obtained by averaging images by groups of 3. 
The final image set obtained this way was much easier to measure than the 
original single frames. In the following discussion of Pic du Midi data we call 
"images" the members of this final set. 

The arc is present on 16 images, and a photometric profile was derived 
for each. Given the poor signal-to-noise ratio (SNR) the curves were later 
average by four, obtaining the result presented in Fig. [5j The profiles have 
been trimmed in order to exclude the region in which the measurement ar- 
eas (main or background patch) are contaminated by the photosphere. In 
other words, only the portion corresponding to the arc projected on the sky 
is plotted. In the case of the Tourelle telescope, the background strongly 
contaminates the signal when the planet is about halfway into the ingress. 
We thus had to be very conservative concerning the fraction of arc to be 
considered. 

Despite the averaging, the result is still affected by a significant noise, 
as indicated by the size of the vertical bars in Fig. |5| representing the 1- 
sigma level. For this reason, small scale fluctuations in the brightness profile 
probably do not have a strong statistic significance if considered separately. 
However, some general trends are present and repeat from one curve to the 
other. For example, a general slope of the curve bottom is present. Also, 
at large values of / (lowest curves) a maximum in intensity is detectable 
close to the South pole of the planet. Interestingly, the maximum is not 
exactly centered on the pole, a feature that repeats in the other data sets 
described further on. The intensity minimum coincides rather well with 
the position angle of the equator of Venus (bar on the left). Even at later 
times, represented by the curves at the highest values, the trend of growing 
brightness toward the pole is present. In all cases, if we neglect the extreme 
tips of the arc very close to the photosphere, the aureole flux is less than 
about 2% the photospheric value, and is detectable down to about 0.2%. 

3.2. Tradate (20 cm Schmidt Cassegrain) data 

L. Comolli observed the event by a digital video-camera, with exposure 
times variable between 1/2000 and 1/50 sec. Only the shortest exposures 
have been used, since they do not present saturation on the solar photosphere, 
thus allowing a normalization equivalent to those presented in the other cases. 
Each analyzed image is the result of the sum of the best 200-300 frames over 
~1 minute of recording. 
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Figure 5: Photometric profile of the aureole during ingress, obtained at the Tourelle 
telescope, Pic du Midi. The plotted value represents the flux coming from a segment of 
1 arcsec of the aureole, normalized to the flux coming from a 1 arcsec^ surface of the 
solar photosphere at 1 Venus radii from the solar limb (measured on the same set of 
images). The horizontal axis shows the position angle along the Venus limb, with the 
usual conventions in the equatorial reference: origin of angles at the celestial N, increasing 
toward E. The two vertical lines represent the position angle of the equator (left, "E" label) 
and the South pole (right, "S"). The flux from the aureole increases with time. Upper 
curves, with the planet about halfway across the solar limb, cover a short angular range 
since their extremes are strongly contaminated by a high background. Labels indicate the 
values of / associated to each curve. 
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Figure 6: Photometric profile of the aureole during ingress, obtained by L. Comolli, nor- 
malized as the other image sets. Vertical bars indicate the entity of the 1-sigma uncertainty 

for each data point. The left bar represents the position angle of the equator, while the 
right bar refers to the South pole. Labels indicate the values of / associated to each curve. 
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Both ingress and egress were imaged and analyzed, using a step of 4.5 
degrees over the Venus hmb. The measured flux distribution during egress is 
well consistent with DOT data (presented below), but limited by a far lower 
resolution. For this reason they are not plotted here. More interesting is the 
comparison of the ingress luminosity profiles to the results obtained by the 
Tourelle Solar Telescope images. 

Five profiles were obtained and are presented in Fig. [6j They confirm, 
from both the point of view of the distribution and the fiux values, the data 
obtained at the Tourelle telescope . Discrepancies around 15% at f~0.5 in the 
faintest part of the aureole appear small when taking into account the wide 
difference in equipment and site, however they clearly hint to the difficulty 
of obtaining clean and appropriate signals at such small solar elongations. 

3.3. Themis Solar Telescope - Pico del Teide Observatory, Tenerife, Spain 

Themis data are complementary to the Pic du Midi observations, since 
from this site it was the exit phase to be imaged, both through a narrow-band 
H-alfa filter and a filter centered at 430 nm. The image sampling is lower 
(Venus radius ~ 70 pixels) but a large image rate is available (1 image each 
1.06 seconds). The sequence was thus divided into different sets of images, 
containing 20 to 60 frames each. The most populated sets correspond to the 
largest values of /, containing an extremely faint aureole slowly fading into 
the background. 

Inside each set, each image was re-sampled with a re- alignment of the disk 
of Venus relatively to an image chosen as reference. An automated procedure 
based on the computation of cross-correlations by Fast Fourier Transforms 
was used, providing the shift to apply to each image. Small relative image 
drifts are thus minimized and the SNR of the faint arc is preserved during 
the summing procedure that yields a final image for each set. 

In this way, 50 images were derived from the sum of the frames inside 
each set, and measured as those of the previous section. The normalization 
procedure adopted was also the same. Due to poor seeing, a width of 7 
pixels was needed to include the entire fiux. On the other hand, the good 
SNR allowed to use a step of 2.25 degrees (corresponding to 2.75 pixels), 
in order to avoid losing information about small scale variations. However, 
the inspection of the resulting curves did not reveal any clear feature at that 
resolution level, much smaller than the arc width smearing due to turbulence. 
A smoothing by running box on three points was thus applied, and the final 
curves are presented in Fig. [7| The overall brightness of the aureole results to 
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Figure 7: Photometric profile of tlie aureole during egress, obtained at the Themis solar 
telescope, Tenerife. The intensity is normalized as in Fig. [5j The two vertical lines 
represent the position angle of the equator (right) and the South pole (left). The curves 
have been split into two families: the upper panel refers to the beginning of the egress, 
while the bottom one, with expanded vertical scale, refers to the central and final phases. 
The labels indicate the values of / associated to the corresponding curves. 
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be 5% of the photospheric reference at maximum, with the shortest measured 
arc. The results at Themis shows some features that are symmetric to those 
observed at the Tourelle telescope during the entry. Here, the brightness 
increases from the equator (right bar in the plot, label "E") toward the pole 
(left bar, "S"). Again, when the polar region is observable projected on the 
sky, a maximum in brightness is detected close to the pole, at about 25 
degrees from it. This is the last portion of the aureole that remains visible. 
Another interesting feature is the relative maximum at position angle ~200°, 
(latitude ~-65°) which is constantly present during most of the evolution of 
the arc profile. 

3.4- DOT telescope - Roque de Los Muchachos observatory, La Palma, Spain 

For several reasons, DOT data constitute the most valuable dataset mea- 
sured for this work, thanks to the optimization of the instrument for solar 
observations ( [Bettonvil et al.[ 2003; Rutten et al. 2004). First of all, the 



image sampling (0.07 arcsec/pixel) provides a comfortable scale for measure- 
ments, resulting in a Venus disk with a radius of 415 pixels. Also, images 
in four narrow bands (usually employed for solar studies) are available, all 
having a FWHM of 1 nm or less (Table [T]). Four cameras were working si- 
multaneously in the four bands, taking every minute an image burst of 100 
frames at 6 frames/sec. For our measurements, all frames of each burst have 
been aligned and summed up to obtain a single image each minute. Due 
to very strong turbulence the alignment process has been very critical, and 
was obtained by applying a Sobel edge enhancement filter with a threshold 



and then using the Hough transform (Yuen et al. , 1989) to find the precise 



position of the disk of Venus in each image. 

Fluxes were obtained for both the aureole and the background with the 
same method as with the previous data set. A step of 2 degrees along the 
Venus limb was used for measurements. We show in Fig. [8] the photomet- 
ric profile referred to the G band. This result is very similar to both the 
blue and red continuum images (Figures [9 10). Only in the Ca II H band 



(Fig. 11) some differences appear under the form of additional peaks in the 
profile. Although a signature of high-contrast chromospheric structures in 
their refracted image cannot be completely ruled out, those features are most 
probably related to the inner corona extending outward from the solar limb 
and contaminating the signal of the aureole, since its flux cannot be easily 
excluded from the measurement areas. 
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Figure 8: Photometric profile of tfie aureole during egress, obtained at the DOT solar 
telescope, La Palma, in the G band. The intensity is normalized as in Fig. [Sj The two 
vertical lines represent the position angle of the equator (right) and the South pole (left). 
The curves have been split into two families: the upper panel refers to the beginning of 
the egress, while the bottom one, with expanded vertical scale, refers to the central and 
final phases. The time interval between each curve is 60 seconds. For an easier visibility 
of the bottom panel, each curve has a step of 7x10"'* added in intensity relatively to its 
lower neighbor. The labels indicate the values of / associated to the corresponding curves. 
Two levels of uncertainty (a) related to a position angle in the range 150-250° are given 
in the two panels. The upper one is applicable to values of / <0.4, while the other one 
refers to /=0.7. 
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Figure 10: As in Fig. |8]in the red continuum at 655 nm. Relative to the other wavelength 
bands, here 3 curves are missing due to corrupted data files. 
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Figure 11: As in Fig. |8]in the Call band at 396.8 nm. 
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However, as far as the main features are concerned, they are very similar 
in all the four bands. For example, in all of them the shape of the profile soon 
after the third contact (first curves in the upper panel. Fig. [s]) has a similar 
shape, showing a bump at high negative latitudes. This aspect is consistent 
with the observations obtained at the Themis telescope. This maximum of 
luminosity at later phases is located between 160 and 210 degrees of position 
angle, and appears to have a double peak. 

Compared to the Themis data, the DOT curves show an overall lower sig- 
nal. Their internal consistency and the difficulties encountered in measuring 
the Themis images affected by poor seeing suggest that the DOT measure- 
ments are probably much more reliable. Also, the low-resolution images by 
ComoUi confirm the flux values obtained for the DOT. For these reasons, this 
data set is chosen for the detailed analysis in the following. 



4. Lightcurves and colors 

In order to study a possible wavelength dependence of the aureole and its 
variation in time, we used the DOT images for the egress phase. The flux was 
measured on three arc segments centered on the pole, on the brightest part 
of the aureole at position angle 195° (position A) and on the faintest part at 



218° (position B) as shown in the scheme of Fig. 12 Point B corresponds to 
a region at latitude 45° S. 

The measurements were made on arc lengths of 10 degrees, in the G band. 



blue and red continuum. The results are shown in Fig. [T3j On all the three 
plots, different wavelengths appear to follow a similar trend. This is expected 
especially for the blue continuum and the G band, which are very close in 
wavelength. 

The fastest fall of flux in red continuum beyond t~12 minutes is a feature 
common to all the light curves, suggesting a certain chromaticity of the 
aureole. However, the related error bars are very large in this interval. At 
the same time, images around t ~ 9-10 min are affected by a momentary 
seeing worsening, so even the color differences at these epochs are to be 
considered with caution. Also, the divergence of the blue and G fluxes beyond 
t ~ 11 min is almost certainly spurious (due to the similar wavelength) and 
further underlines the difficulty of comparing very low fluxes. 

In conclusion, we cannot flrmly state that the images at our disposal 
reveal a signiflcant departure of the aureole color from the solar spectrum. 
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Figure 12: Geometry at / = 0.5 during the final phases of the event. The large dots 
indicate the points that where chosen for studying the variation of the arc brightness in 
time. A corresponds to the brightest portion, and B to the faintest. 

5. Modeling the aureole 

5.1. The refraction model 

The aureole observed during Venus transits can be explained by the re- 
fraction of solar rays through the planet outer atmospheric layers. The rays 
that pass closer to the planet center are more deviated by refraction than 
those passing further out. The image of a given solar surface element is 
flattened perpendicularly to Venus' limb by this differential deviation, while 
conserving the intensity of the rays, i.e. the brightness of the surface element 
per unit surface. This holds as long as the atmosphere is transparent, i.e. 
above absorbing clouds or aerosol layers. 

The refractive deviation of light is related to the physical structure of 
the planet atmosphere. The formalism of this deviation has been studied by 
Baum and Code (1953), for the cases of stellar occultations by planets. Their 
approach assumes that (i) the local density scale height H of the atmosphere 
is constant and much smaller than the planet radius R^, (ii) the atmosphere 
is transparent, {Hi) and it has spherical symmetry. 

This approach is still valid in our case, but has to be modified to ac- 
count for the finite distance and size of the Sun. In the following we will 
consider that the physical properties of the transparent atmospheric layers 
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Figure 13: Relative intensity of the arc brightness in the DOT images, normahzed as in 
the previous plots, as a function of time, where the origin of time corresponds to the 
third contact (at 11:07:25 UT). The three panels refer to the Pole and to points A, B 
(in the order, from top to bottom). The G band, red continuum, and blue continuum 
are represented by the solid, dotted and dashed line respectively. The model fit is the 
thickest solid line in the top and bottom panel. At the B point, H=3.1 km and Ar=26 km 
provide the best fit, while at the pole the curve corresponds to H=4.8 km and Ar=38.5 km 
(see Sect, [sjfor details). The cross at t~15 min represents an isolated red measurement 
available. Error bars are very similar for all the three bands. 
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vary smoothly with altitude, and remain about constant over a scale height. 
Note that the assumption <^ i?^ is valid here, as H is smaller than the 
atmospheric layer of Venus, in turn -C the planet radius. Note also that 
the locally spherical symmetry is achieved for Venus' atmosphere. As to the 
transparency assumption, it must be dropped when the rays are going deep 
enough for the atmosphere to become opaque. 
We model the aureole brightness as follows: 

(1) We consider a surface element dS on the Sun. A ray emitted by 
dS will reach the observer O after being refracted by an angle u: in Venus's 



atmosphere (see Fig. 14). Note that we adopt here the convention w < 0. 



To each element dS corresponds an image dS' that will appear as an aureole 



near Venus's limb (Fig. 15) 



During its travel to Earth, the ray passes at a closest distance of r from the 
planet center. Furthermore, dS projects itself at a distance from Venus' 
center C, as seen from O. Note that is an algebraic quantity which is 
negative if dS and dS' project themselves on opposite side relative to the 



planet center C, and positive otherwise, see Figs. 14 15 



(2) The refraction angle oj is given by Baum and Code (1953): 

-v{r)^j2T^r/H, (1) 



where v is the gas refractivity, decreasing exponentially with r. This quantity 
is related to the gas number density nhy u = K ■ n, where K is the specific 
refractivity. For CO2, we hav^i^ = 1.67 x 10^^^ m^ molecule"^. 

(3) It can be shown that surface element dS' is radially shrunk with 
respect to dS by a factor $ = 1/ [1 + D ■ {duj /dr)], where D is the distance 
from Venus to Earth. Note that since the atmosphere is assumed to have a 
constant density scale height, doj/dr = —u/H, thus: 

$ = ^ (2) 

(4) It is convenient to take as a reference radius the closest approach 
distance ri/2 corresponding to $ = 0.5 (so so-called "half-light radius" in the 
stellar occultation context). Thus UJ112 = —H/D. Re-arranging the various 



equations given above, one finds K ■ 721/2 = \ I (27rri/2-D2) , where 



n 



1/2 



see http://www.kayelaby.npl.co.uk/generaLphysics/2_5/2_5_7.html 
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Figure 14: Geometry of the refraction of solar rays by Venus' atmosphere, where D' 
(resp. D) is the distance of Venus to the Sun (resp. Earth). A ray emitted from the solar 
surface at dS, is deviated by an angle uj (here negative) before reaching the observer O, 
who observes the images dS' just above Venus' limb (the aureole). The observer sees dS 
projected at algebraic position in the plane going through Venus and perpendicular to 
the line of sight. The origin of is at Venus' center C and increases upward, see also 



Fig. 15 All sizes and angles have been greatly increases for better viewing. 
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Figure 15: Venus (dark gray disk) observed from Earth, partly against the solar disk and 
partly against the sky (black background). Each solar surface element dS surrounding 
ri has a refracted image dS' of length I and width dr', caused by Venus atmospheric 
refraction. The image dS' has the same surface brightness as dS if the atmosphere is 
transparent. See text for details. 



25 



is the gas number density at ri/2- The latter equation permits to derive the 
numerical value of the half-light radius, once an atmospheric model of the 
planet is given, that is, once the density profile n{r) is specified. 
(5) From Eq. [l| we derive: 



H 
D 



' 1/2 



(3) 



Simple geometrical consideration show that we also have: u = —k (r — rj) /D, 
where k = 1 + D/D' [D' being the distance of Venus to the Sun). Combining 
this equation with Eqs. [2]j3l we obtain: 



This is the Baum and Code formula (apart for the correcting factor k, 
which is equal to unity in the original formula since D' = +oo for stellar 
occultations). 

Thus, for each value of rj, we can calculate the corresponding value $(rj), 
using a classical Newton iterative numerical scheme. 

For each value of rj, we can also calculate the corresponding closest ap- 
proach radius r by combining Eqs. |2] and [3] 

r = n/,-H. \og(^-l] (5) 



To calculate the flux dF received from an aureole element of surface dS' 



with length / and width dr' (Fig. 15), it is enough to note that the surface 



brightness of dS and dS' are the same if the atmosphere is transparent, i.e. 
that dF = SQ{ri)ldr', where SQ{ri) is the flux received from a unit surface 
on the Sun at rj (taking into account the limb-darkening effect). Thus, we 
can re-write this equation as dF = SQ{ri)^{ri)lidri. 

The aureole is not radially resolved, so we have only access to the flux 
integrated along r^, i.e. to: 

' ^'"^^^ SQ{nmn)-k-dn (6) 



2, mm 



The lower bound of the integral, rj^min, corresponds to the value of r 
corresponding to an opaque cloud or aerosol layer at altitude icut- The upper 



26 



bound of the integral, max, corresponds the solar limb, beyond which no 
more photons are emitted. 

By applying this model, in the following we will determine the scale height 
H and the half occultation radius relative to slanted opacity r ~1 (Ar = 

— best reproducing the observations. In general, different portions of 
the arc can yield different values of these parameters, thus providing a useful 
insight of variations in the physical properties of the Cytherean atmosphere 
function of latitude. 

6. Derivation of the physical parameters 

We used the equation (|4]) for modeling the refraction in the atmosphere 
of Venus and the flux in the observed aureole. Essentially, the mathematical 
model depends upon two parameters: the altitude of the half-occultation 
level relative to an underlying totally opaque layer (Ar) and the optical scale 
height of the atmosphere H, second free parameter of Eq. |4j 
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Figure 16: Relative intensity of the arc brightness for H=4.8 km, Ar^38.5 km as derived 
from the model. Labels correspond to the / fraction of disk overlap. 

Since our brightness measurements are referred to different epochs, the 
model must fit at the same time not only single-epoch profiles, but also their 
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evolution. Although free parameters of Eq. |4| ri/2 and H are quantities that 
are not physically independent since they are both related to local physical 
structure of the mesosphere of Venus. We will thus check afterwards the 
physical consistency of our results. 

The source of refracted rays (the solar photosphere) is computed as a 



smooth function of the radial photometric profile of the Sun (given by Hes- 



troffer et al. 1998), whose parameters are determined by a fit to the profile 



measured on the solar photosphere imaged close to the Cytherean disk. 

Fig. 16 shows a typical outcome of the model for different phases / of 
Venus ingress or egress, considering the atmospheric properties to be constant 
all along the planet limb. When Venus appears largely overlapped to the solar 
disk, the refracted flux is dominated by the arc extremes, where it approaches 
the solar limb. This is a consequence of a very small bending of the light path 
occurring high in the atmospheric layer involved. When f > 0.5, the regions 
of the photosphere contributing to the arc brightness are all farther away 
than the Venus disk center - so the larger deviation associated to a deeper 
pass of the solar hght inside the atmosphere is needed. In such situation, 
the arc extremes are fainter than the center, since their brightness sources - 
being on the opposite side of the planet center - are on the strongly darkened 
photospheric limb. For higher values of / the signal rapidly fades. The cut- 
off at zero flux is a feature of the model due to the opaque layer, blocking 
the refracted rays that cross the atmosphere below a given altitude. 

Diminishing H determines a stronger refraction that increases the amount 
of light coming from regions that are farther away from Venus on the solar 
photosphere. The slope of the flux decrease that is observed when Venus 
exits the photosphere (i.e. when / is increasing) will also be less steep (Fig. 



17|). 

Another feature of the model outcome is the symmetry relative to the 
axis joining the center of the disks of Venus and the Sun. It is thus clear that 
the asymmetries and features visible in the measurements must be related 
to latitude variations of the physical properties in the refracting atmospheric 
layers. 

We searched for the values of Ar and H that better fit the time variations 
reproduced in Fig, 13 , separately for the Pole (or, equivalently, the point A) 
and the point B. The procedure we adopted started with a search for the 
value of Ar, which mainly affects the amount of aureole flux at small / (i.e. 
when Venus is nearly completely on the solar disk). Then, H is varied in 
order to obtain the right slope for the flux variation in time. The procedure 
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is repeated until a convergence to the observed data is obtained. 

The central portion of the arc, for which the longest time series is avail- 



able, is best fitted by Heg=3.1 km and ATeq=26 km (Fig. 13). This region 
corresponds to intermediate latitudes (~45° - point B). 

We then fitted the Pole area, obtaining Hpo;e=4.8 km and Arpo/e=38.5 km. 
This determination appears less satisfactory at low / values (beginning of 
egress). If a better fitting of this part of the curve is forced (in particular by 
increasing Ar) the central plateau of the curve becomes more extended, and 
the fit is worsened in the rest of the domain. We thus preferred to privilege 
the overall better fit. This choice appears to be reasonable when compared 
to the expected altitude difference of cloud top and upper haze between the 
equator and poles (see next section). 

In fact, if refraction occurs at similar absolute altitudes in the atmosphere, 
variations of Ar should only correspond to changes in the altitude of the 
opaque layer, i.e. a slanted opacity (optical thickness) r ~1. Our best fit 
values indicate that the higher flux received at the polar latitudes on Venus 
can be related to a significantly lower aerosol extinction towards the poles, 
thus allowing more refracted rays to reach the observer. This contribution is 
blocked at the equatorial latitudes by the higher clouds. 

An indication of the formal accuracy of the determination of Ar and H 
can be obtained by studying the fit sensitivity to changes in the parameters. 
We have verified that a change of ~3 km in Ar and ~100 m in H are sufficient 
to displace the fitting curve of an amount exceeding the extension of the error 
bars all along the lightcurves. The rapid variations in the aureole behavior 



following changes in the parameters are illustrated in FigjT7| In particular, 
our results appear extremely sensitive on H. As a further consistency check, 
we can estimate Ar from Eq. [3] by using a value of H~4 km and assuming 
the maximum deviation of light rays (corresponding to a highest / at which 
the signal from the aureole is detected) to be of the order of u=40 arcsec. In 
this case Ar~35 km, very close to the results given above obtained by a full 
model fit to the observed data. 

Eventually we can estimate the visibility of the aureole in terms of visual 
magnitude. We use the simple law for the profile of the solar limb darkening 



(normalized to the Sun's center brightness) by Hestroffer et al. (1998): 



/(/i) = 1 - n (1 - /i") (7) 

in which yU = (1 — r^)°'^. Here, r represents the distance from the center of 
the Sun, normalized to the solar radius. At a wavelength in the V-band range 
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Figure 17: EflFects of the variation of the physical parameters controlling the aureole flux. 
In the top panel, H(.y=4.8 km is kept constant, while the three simulated light curves 
correspond to different values of Argg (indicated by the labels). On the other hand, in the 
bottom panel Areg=38 km is constant and H is varied. 
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(A=579.9 nm) this law fits the observed profiles when u=0.85 and a=0.8. In 
the case of our reference area the value of r is r„=0.969, so It,=0.428. In 
this region so close to the solar limb the profile approximation can be less 



accurate, but it should not be erroneous by more than 2% (Hestroffer et al 



1998D . 

From the computed I^ and using the magnitude of the Sun (V = -26.71) 
we obtain a surface magnitude per arcsec^ for the reference element m^ = 
— 10.87. By observing the plots, we can assume that the aureole has been 
easily seen and imaged when its brightness was ~10~^ times the surface 
brightness of the reference area (1 arcsec^ at 1 Venus radius from the solar 
limb), i.e. 5 magnitudes fainter. We can thus consider that the typical 
magnitude of an arc of aureole 1 arcsec in length was -5.9. 

One should note, however, that ample variations around this value are 
present, both along the arc and over time. 



7. Comparison to Venus Express observations 

Venus' mesosphere, which extends typically from 70 to 110 km, appears 
to be a transition region between the troposphere (0-70 km), dominated 
by the ~4 Earth days retrograde zonal superrotation, and the thermosphere 
(110- 250 km) in which solar EUV and diurnal temperature contrasts drive a 
predominantly subsolar-to-antisolar (SS-AS) circulation. Monitoring of the 
mesopheric thermal structure has revealed that this region is characterized 
by a strong temporal variability, whose origin remains poorly understood 



(Bougher et al. 


1997 


2006 


Lellouch et al. 




1997 


2008 


Sandor and Clancy 


2005 


Gurwell et al. 




2007 




Widemann et al. 


2008 


). Upward from cloud top 



near 70 km, zonal winds decrease with height while thermospheric SSAS 
winds increase. The region further serves as the primary photochemical re- 
gion of the Venus atmosphere and heterogeneous chemistry revealed by haze 
and cloud particles production. Although present in small amounts, trace 
gases are involved in important complex chemical cycles. Photochemical 
reactions between CO2, SO2, H2O and chlorine compounds lead to the for- 
mation of sulfuric acid, which is the main component of the cloud and upper 



haze, associated with low H2SO4 vapor pressure (Zhang et al. , 2010). 



7.1. Cloud and haze structure 

At the limit between the troposphere and the mesosphere, the upper 
cloud structure reveals distinct dynamical structure in the latitudinal direc- 
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tion, with convective, wave-dominated zonal flow in the lower latitude range 
near Venus's equator, and a significant transition to less convective banded 
flow between 45° up to about 70°-80°. The temperature field forms a bulge 
of cold air at 60°-80° latitude called the "cold collar region", which verti- 
cally extends up to 75 km. The torus-like structure encloses a vast vortex 
several thousand kilometers across, with a slower rotation period (~2.5 to 
~2.8 Earth days) and significantly depressed cloud top altitude, at about 



62-65 km instead of 70-74 km closer to the equator (Piccioni et al. , 2007). 



The structure of the cloud tops is especially poorly investigated since it falls 
between the altitude ranges sounded by solar/ stellar occultations and that 
studied by descent probes. Recent analysis has been performed on board 



Venus- Express (Svedhem et al. , 2007) based on depth of CO2 bands at 1.6 
yum measured by VIRTIS (Ignatiev et al. 2009), and CO-I-CO2 gaseous ab- 
sorption in the 4.5-5 fim range using VIRTIS and VeRa temperature profiles 



(Lee et al. 2010) while Luz et al. (2011) brought the first extensive charac 



terization of the vortex dynamics and its precession motion. Although the 
absolute cloud top altitude poleward and equatorward of the polar collar 
differs when applying the two modeling techniques (63-69 km / 74 ± 1 km 
in Ignatiev et al., 2009 ; 62-64 km / 66 km in Lee et al., 2010), their relative 
difference is comparable (1-2 scale heights). The latitudinal extent of this 
polar depression matches well the brightest portion of the aureole observed 
southward of 65°-70° near the South polar limb in DOT telescope data (see 
Fig.Sa of Ignatiev et al. 2009). Interestingly, similar agreement can be traced 
back to Russell's drawings of 1874 (See Link, 1969 and Fig. 1 of Pasachoff 
et al, 2011). 

7.2. Haze aerosols 

The Venus upper haze (70-90 km) was first evidenced by measurements 
from Pioneer Venus orbiter limb scans at 365 and 690 nm at northern mid- 
latitude. It is mainly composed of submicron sulfuric acid (H2SO4) aerosol 
particles with typical radii from 0.1 to 0.3 fim (Lane and Opstbaum, 1983 
; Sato et al., 1996). Wilquet et al. (2008) also present the first evidence 



for a bimodal particles distribution, similar to the two modes in the upper 
clouds, at the latitudes probed during the observation of solar and stellar 
occultations by ESA's Venus Express, with typical radii of mode-2 particles 
between ~0.4 and lyum. These measurements were performed close to the 



polar regions. From the results of Wilquet et al. (2008) the absorption at dif- 



ferent wavelengths as a function of the altitude can be deduced. In general. 
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the absorption levels are found to be 6 ± 1 km higher in the visible domain 
(SPICAV-IR at 757 nm) than at 3 /xm. 

As an integrated aerosol optical depth ~1 is measured a few density 
scale heights above the cloud tops, we note that the slanted geometry of 
the aureole must reach half occultation level ri/2 well above the upper haze. 
Quantitatively, the altitude of the aureole's half-occultation level in the polar 
region will be found by adding the value of Arpoie=38.5 km to the altitude 



where aerosols slanted opacity r ~1. Recent VEx/SOIR results (Wilquet et 



al. , 2011, personal communication) place that altitude at 73 ± 2 km in the 



3/xm band, to be further increased by 6±1 km to retrive its value in the visible 
domain, as above. The final sum thus yields ri/2 ~ 117.5 ±4 km at the pole. 
Most recent VEx/SOIR results integrated over four years of observations 
(Wilquet et al., 2011) further indicate an altitude of longitudinally averaged, 
integrated aerosol optical depth ~1 increasing toward the equator, with an 
altitude of 81 ± 2 km for latitudes between 35°S and 55°S. Therefore, as 
Aieq =Ti/2—Tcut IS significantly smaller than Ar^o/e the resulting value of 
Ti/2 at mid- latitude differ from our calculation in the polar region by about 
two scale heights (Table 2). We can thus consider that our results only 
show marginal latitudinal variation of the altitude for the half-occultation 
level along the terminator, although in a region of important temperature 
variation which need to be independently assessed. The same SPICAV/SOIR 
results show that, under reasonable assumptions, the refraction index at 
visible wavelengths is fairly constant, thus confirming that the process is not 
a source of relevant chromatic effects in the aureole. 



7.3. Temperature structure 

Only scarce measurements of vertical temperature profiles have been per- 
formed above 100 km altitude, where temperature fields, especially in the 
polar region, as well as their time variability, are still debated (see, e.g., 
Vandaele et al. 2008, Clancy et al. 2008, 2011; PicciaUi et al. 2008, 2011). 
Inverted equator-to-pole temperature gradient on isobaric surfaces above 75 
km, were first reported by NASA's Pioneer Venus Infrared Radiometer and 
radio occultation measurements (Taylor et al., 1983; Newman et al., 1984), 
Venera-15 Fourier Spectrometry data (Zasova et al., 2007) and more recently 
ESA's Venus Express, e.g. Grassi et al. (2008) ; Piccialli et al. (2008) us- 
ing VIRTIS -M observations ; Tellmann et al. (2009) based on VeRa radio 
occultations on board VEx. The collar region also divides the atmosphere 
vertically. Below the collar, the atmosphere cools with increasing latitude. 
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Table 2: Results of modelling and retrieved half-light altitude ri/2, scale height H (in km) 
and temperature along the Lomonosov's arc, assuming a locally isothermal atmosphere. 
Explanations of symbols in the text. Tjso ~ m(z) g(z) H/k. Tmod (K) from Hedin et al. 
1983. 



Point 


Lat. 


Ar 


H 


^cut 


ri/2 


Tiso (K) 


Traod (K) 


A 


68°S 


38.5±3 


4.8±0.1 


76±2 


120.5±4 


212 


138-190 


B 


45° S 


26.0±3 


3.1±0.1 


81±2 


113.0±4 


137 


143-185 



Above, the temperature gradient is reversed, as diabatic heating or dynam- 
ically controlled processes could be responsible for the observed structures 
(Schubert et al., 1980, Tellmann et al., 2009). At the altitude of the collar 
and closer to the poles, the atmosphere is almost isothermal and also much 
warmer than the surrounding areas. Vertical temperature profiles are fairly 
constant with altitude between around 50 and 100 km, above which they 



begin to rise sharply in the mesosphere, see e.g. Fig. 2 in Mueller- Wodarg 



et aL| ( |2006D . 

Retrieving the temperature using the density scale height should be made 
with caution. The transition between stable and adiabatic regions is fairly 



smooth at equatorial latitudes, but abrupt at middle and polar latitudes (Tell- 



mann et al. , 2009). A significant vertical temperature gradient is a potential 
issue when trying to compare the density scale height to the temperature 
scale height (see the discussion related to Titan's mesosphere refractivity 



m 



Sicardy et al. , 2006). However, we consider the mid-occultation level 
probed by the arc in the ~115 km region (see Table 2) as isothermal. This 
is an acceptable hypothesis since the diurnal average of dT/dz on the back- 
ground atmosphere (Hedin et al., 1983) is negligible, although significant 
local-time and daily variability has been retrieved (Clancy et al. 2008, 2011; 
Sonnabend et al. 2008, 2011). In particular, the elevated value of Tj^o = 
212K at 120. 5±4 km at 68°S is in agreement with high kinetic temperatures 
derived from heterodyne mid-infrared spectrum of non-LTE emission at ter- 
minator, see in particular Table 1 of Sonnabend et al. 2011. This result 
suggests a warmer or rapidly variable temperature condition in the altitude 
range. 

7.4- Physical interpretation 

In order to check the consistency of our results with the known physics 
of Venus atmosphere, we can compute the expected value of the density 
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scale height H at the altitude of the refracting layer. The temperature scale 
height is defined a.s H = kT{mg)~^, where k is the Boltzmann's constant, 
while T, m and g respectively represent temperature, molecular mass and 
gravitational acceleration at the altitude considered. As discussed above, we 
consider a locally isothermal atmosphere. For our computation we derive the 
value of g as a function of altitude. We derived an altitude for the opaque 
layer from SOIR aerosol transmittance measurements near 3 /xm (Wilquet 
et al., 2011 ; Wilquet, personal communication), and added 6 ± 1 km for 



transferring to the visible domain. We then adopt the model by Hedin et al. 



(1983) describing the atmospheric structure on Venus at noon and midnight, 
around the equator. We consider that an appropriate temperature estimation 
can be obtained by averaging the two profiles. Differences may be present, 
since we observe directly at the terminator, but the temperature data by 



SOIR obtained at a similar geometry (Mahieux et al. 2010) do not show 



significant discrepancies relative to the model. CO2 is the dominant gas on 
the dayside below 160 km and on the nightside below 140 km, replaced at 
higher altitudes by O. Mean molecular mass remains close to 44 a.m.u. below 



130 km, where CO2 dominates, and decreases towards higher altitudes (Hedin 



et al. , 1983 Mueller- Wo darg et al. , 2006). Although day-night differences 



in composition are substantial, we also may consider longitudinally averaged 
values at the altitude considered. Concerning m, since the layer is in the 
heterosphere, a decrease due to fractionation has to be taken into account, 
due to the decrease in CO2 relative to other components, in particular atomic 
oxygen O, molecular nitrogen N2 and He. As a result, at 120 km we use an 
average molecular mass of 42.9, and a day-night longitudinal average T=164 
K, obtaining H=3.7 km. An equal value is obtained at 115 km. Both are 
in remarkable agreement with the average one obtained from the aureole 
measurements, H=3.95 km (see Table 2). Finally, we noted that our results 
only show marginal latitudinal variation of 1-2 scale heights for the altitude 
for the half-occultation level along the terminator, near 117 km. 



8. Conclusions 

Our results represent the first successful model of the aureole of Venus, 
observed during solar transits. We are able to reproduce the main features 
of the lower mesosphere observational constraints : longitudinally averaged 
thermal structure near 117 km, slanted opacity of aerosols and their merid- 
ional variation, with consistent physical parameters. 
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The results obtained by this first set of measurements are encouraging 
and suggest that more accurate planning will produce more precise data. 
Improvements are also possible concerning the refraction model, which could 
be completed by a more realistic, gradual transition from transparency to 
absorption at the cloud deck level based on most recent results on the latitu- 
dinal distribution of upper haze aerosols. Also, the contribution of scattered 
light in addition to refracted light could be included. In fact, the absorp- 
tion curves in Wilquet et al. (2008) Fig. 3 are wavelength dependent, i.e. 



the opaque layer altitude in our model can be function of the color. As a 
result the aureole should not be perfectly "grey", suggesting that future ob- 
servations capable of accurate multicolor photometry should investigate this 
possibihty.Even a grey absorption would cause the solar flux in the aureole 
to look dimmer than without aerosols, which would be wrongly interpreted 
as a drop in scale height, i.e. in local temperature at tangent point altitude, 
in the inversion code. 

In the two hemispheres, the cold collar temperature structure and the 
polar regions are very similar (Tellmann et al. , 2009), so the latitudinal vari- 
ations obtained for the 2004 transit around the South polar region might 
tentatively apply to the North polar data we expect to acquire in June 2012. 
Recent interpretations also attest a general N-S symmetry in the latitudinal 
variation of aerosols extinction (Wilquet et al. , 2011). Given the difference 
in brightness between the arc and the solar photosphere, accurate measure- 
ments of this phenomena are challenging. The degree of turbulence, the 
magnification, the amount of light scattering in the optics are all elements 
that contribute in determining the effective visibility of the aureole in a given 
instrument, and the accuracy of the measurements. The transit of June 2012 
will be our last opportunity for observing the aureole using this spatially 
resolved technique, until the next pair of transits of Venus, which will be in 
the ascending node, on 11 December 2117 and 8 December 2125. 
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